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Abstract 

The orbital angular momentum is one of the least understood of the spin 
characteristics of a proton. There are no direct ways to model L z . How- 
ever, the J z = \ sum rule includes an angular momentum component and 
can provide indirect access to properties of L z . One of the other unknowns 
in the sum rule is the gluon polarization, AG. We can define the gluon spin 
asymmetry in a proton as A(x,Q 2 ) = ^T^pv • This can be written as a 
sum of a Q 2 -invariant piece, A (x) and a small C} 2 -dependent term, e(x, Q 2 ). 
The x-dependence of Aq can be calculated and a suitable parametrization for 
e(x,Q 2 ) can be made to estimate this asymmetry. When combined with the 
measured unpolarized gluon density, G(x, Q 2 ), this provides a model indepen- 
dent prediction for AG(x,Q 2 ). This eliminates one unknown in the J z = \ 
sum rule and allows a reasonable estimate for the size and evolution of the 
orbital angular momentum of the constituents, L z . 

1 Introduction 



High energy spin physics studies have evolved into three primary areas. The first 
involves helicity components of the constituents of the proton spin. From the initial 
EMC experiments that defined the "proton spin problem", these studies have fo- 
cussed on determining the proportion of spin carried by the valence and sea quarks 
in addition to the gluons. The remainder of the nucleon spin is attributed to the 
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orbital angular momentum. The second important area of spin studies involves 
transversity measurements. These not only involve the transverse motion of the 
nucleon constituents, but the dynamics of the interactions among them. The Sivers, 
Collins and Boer-Mulders functions all contain information about transversely po- 
larized quarks in various polarized states of nucleons. In all cases, this information 
is related to the orbital motion carried by these constituents. Finally, our studies of 
generalized parton distributions constitute an attempt to put these and other spin 
phenomena into a general formalism. In all of these areas, the nature of the orbital 
angular momentum of nucleon constituents plays an important role. 

2 Determination of AG(x) with the Gluon Asym- 
metry 

The J z — \ sum rule can be used to access the orbital angular momentum phe- 
nomenologically. This involves the integrated parton densities. 

J z = ^AE + AG+(L z ) q+G , (1) 

where AE is the total spin carried by all quarks, AG the spin carried by gluons and 
(L z ) q+ c is the orbital angular momenta of the quarks and gluons. Numerous DIS 
experiments have narrowed the quark spin contribution (AS) to within a reasonable 
degree. However, the gluon and orbital angular momentum components of the 
nucleon spin are virtually unknown. 

Experiments are presently underway in various kinematic regions [1, 2] to deter- 
mine this distribution. Others have been proposed to expand these measurements 
to other kinematic regions. [3] Meanwhile, there have been many models assumed 
for AC Our calculation of the gluon asymmetry does not presume any specific 
model for the polarized gluons, but relies on simple theoretical assumptions and the 
measurements of the unpolarized gluon distribution. 

We define the gluon polarization asymmetry as 

A(x,t) = AG(x,t)/G(x,t), (2) 

where the evolution variable t is defined as t = hi[a s {Ql) / a s {Q 2 )]. Since there 
are no overwhelming theoretical arguments favoring any single model for AG, we 
consider a more direct argument for its shape in terms of this asymmetry. The 
conclusions we draw follow from the observation that, in the absence of a "valence" 
gluon, both G(x,t) and AG(x,t) exhibit scaling violations which can be associated 
with measurements resolving radiative diagrams. The diagrams leading to positive 
and negative helicity gluons are the same. This implies that the relative probability 
of measuring a gluon of either helicity does not depend upon t. Thus, to a first 
approximation, the gluon polarization asymmetry is predicted to be at most, mildly 
scale dependent. 
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To construct a theoretical model of the asymmetry, we assume that it has a scale 
independent part, A (x) plus a small piece that vanishes at some large scale. Thus, 
we can write the asymmetry as 

A(x, t) = A (x) + e(x, t) = AG/G. (3) 

Then, AG can be written in terms of the calculated asymmetry A (x) and a dif- 
ference term as AG = A ■ G + e • G. The scale invariant A is calculable and is 
independent of theoretical models of AG. The second term is interpreted as the 
difference between the measured polarized gluon distribution and that predicted by 
the calculated A combined with measurement of the unpolarized gluon density. It 
is reasonable to choose t — to coincide with a typical hadronic scale, Q 2 = ml. 
The calculable part can be found by taking the t-derivative of A (x): 

dt 1 dt dt dt i dt v 1 

We assume at some scale that the quantity e • G is t-independent. Thus, 

dAG 

= (5) 

dt 

Both and ^ are calculated using the DGLAP evolution equations. So 



dAG 

(6) 



A <H " AP Gq ® Aq + AP GG <s AG 



dG 
dt 



P Gq ®q + Pgg ® G 



Since AG has not been measured, equation 6 can be converted into an iterative 
equation for A(x) by inserting AG(x, t) = A(x) ■ G(x, t) from equation 3 into the 
convolution, 

~AP Gq ®Aq + AP GG ® (A n ■ G) 



i-n+l — 



(7) 



P Gq ® q + P GG ® G 

This equation then can be solved iteratively. With a suitable assumption for e(x, t = 
0), the asymmetry A(x, t) can be determined at all x and t values. From the counting 
rules, we bound e(x, t) by e(x, t) < c(t) ■ x{l — x) and require e(x, t) to be decreasing 
at some scale, since its evolution is opposite that of the gluon. Then, the form for 
e(x, t = 0) = x{l — x) n , where n is the power of (1 — x) in G(x) and we assume that 
c(0) = l. 

For the starting distributions in equation 6 and the iterations of equation 7, 
we use the polarized quark distributions outlined by GGR [6] using the CTEQ4M 
and MRST unpolarized distributions. [7, 10] The evolution was performed in LO, 
since the NLO contributions to the splitting kernels, calculated in ref. [8], are most 
dominant at small-x, where the asymmetry is the smallest. Work is in progress to 
ensure that the effects of NLO are not significant for the ratio ^p. The iteration is 
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relatively stable and converges within a few cycles. Some models of G(x) converge 
more uniformly than others. The resulting range of shapes for A(x) generated by 
these distributions is shown in Figure 1. The corresponding extremes of xAG are 
shown in Figure 2. 

The HERMES experimental group at DESY has measured the longitudinal cross 
section asymmetry A\\ in high-p^ hadronic photoproduction. [1] From this and 
known values of — from DIS, a value for A G (x G ) are extracted. Here, xg = s/2Mv 
is the nucleon momentum fraction carried by the gluon. The COMPASS group has 
also measured this asymmetry at a slightly smaller average value of x G . [11] The 
results of these measurements are 

• HERMES: A G = 0.41 ± 0.18 (stat.) ±0.03 (syst.) at < x G >= 0.17 

• COMPASS: Aq = 0.06 ±0.31 (stat.) ±0.06 (syst.) at < x G >= 0.095. 

Our range of asymmetries falls within these values. 

The corresponding calculation of L z and its evolution involves using the J z — \ 
sum rule and the DGLAP evolution equations. 

L z = \~ AS/2 - (A + e)-G. (8) 

From its derivative with respect to t and the evolution equations, the evolution of 
L z is 

^ = -[AP qq ® Ag + AP qG <g> ((A ± e) • G)\/2 - (A ± e)[P Gq ®q + P GG ®G}. (9) 



3 Results and tests for L z and its evolution 

A plot of the L z (x) is shown in Figure 3. Differences in the MRS and CTEQ based 
distributions can be seen from the range shown. The evolved L z to 100 GeV 2 for the 
CTEQ based model is also shown in the figure as a dotted line. The range of results 
for the integral of L z is consistent with those outlined in reference [12] and the plots 
of L z {x) are comparable to those in reference [13]. The evolved L z to Q 2 = 100 
GeV 2 shows a positive component of L z at moderate to large x. This can be tested 
in future experiments. 

All of the results presented here can be tested by three separate experiments 
at DESY (HERMES), BNL (RHIC-STAR and PHENIX) and CERN (COMPASS). 
First, AG can be measured via prompt photon production or jet production. These 
processes yield the largest asymmetries for the size of AG. [2, 9] The kinematic 
regions of STAR and PHENIX can determine A over a suitable range of XBj to test 
this model of the gluon asymmetry. Coupled with additional direct measurements 
of A(x G ) from HERMES, an appropriate cross check of AG(x) and G(x) can be 
made. Since the range of values for ^ encompass both a comparable and larger 
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Figure 3: Range of L z versus x corresponding to the range of asymmetries in Figure 
1. The dotted line is L z evolved to Q 2 = 100 GeV 2 . 

polarized glue than the GGRA model used in ref.[9], all of the asymmetries for 
direct-7 and jet production should be able to narrow the values due to possible 
enhanced asymmetries. 

The COMPASS group at CERN [3] plans to extract A from the photon nu- 
cleon asymmetry, A™ n (xg) in open charm muo- production, which is dominated by 
the photon- gluon fusion process. This experiment should be able to cover a wide 
kinematic range of xq as a further check of this model. The combination of these 
experiments will be a good test of the assumptions of our gluon asymmetry model 
and a consistency check on our knowledge of the gluon distribution in the nucleon 
and its polarization. The corresponding predictions for the orbital angular momen- 
tum L z and its evolution can be measured via deeply-virtual Compton Scattering 
(DVCS) in the HERMES or COMPASS experiments. 

We have outlined a possible method, based upon simple theoretical assumptions, 
to estimate AG(x) and thus provide some insight as to the x and Q 2 dependence of 
the orbital angular momentum. All aspect of this model can be readily tested in the 
experiments of the HERMES, RHIC (STAR and PHENIX) and COMPASS groups. 
The results will bring us closer to understanding the nature of spin in fundamental 
particles. 
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